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Thermodynamics of the exactly solvable spin-electron tetrahedral chain 
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An exactly solvable spin-electron tetrahedral chain, where the Ising spins localized at nodal lat¬ 
tice sites regularly alternate with three equivalent lattice sites available for one mobile electron is 
considered. The system with ferromagnetic interaction between the Ising spins and electrons ex¬ 
hibits an enhanced magnetocaloric effect (MCE) in the limit H/\J\ —>■ 0 when the entropy is very 
small, whereas the system with antiferromagnetic interaction between the Ising spins and electrons 
exhibits an enhanced MCE around the field H/J — 1 when the entropy is sufficiently close to the 
value S/2N = ln[(l + %/5)/2]. We study the thermodynamics of the system in these field regions. 
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I. INTRODUCTION 

Exactly solvable ID models provide an excellent play¬ 
ground for theoretical studies of collective phenomena. 
One such model is a coupled tetrahedral chain, which has 
been exactly solved for the spinless fermion model [l| and 
the mixed Ising-Heisenberg model 0. In this paper, we 
consider another exactly solvable model with tetrahedral- 
chain structure, where the Ising spins localized at nodal 
lattice sites alternate with three equivalent lattice sites 
available for mobile electrons. Note that the particular 
case of this hybrid model with one mobile electron delo¬ 
calized over each triangular plaquette has been recently 
solved to study the ground state, the magnetization pro¬ 
cess and the magnetocaloric effect of the system Q . The 
aim of the present work is to shed light on the thermo¬ 
dynamics of the model in those field regions, where an 
enhanced MCE has been found. 


II. SPIN-ELECTRON TETRAHEDRAL CHAIN 

We consider a tetrahedral chain in an external mag¬ 
netic field, in which one Ising spin localized at the nodal 
lattice site regularly alternates with three equivalent lat¬ 
tice sites available for one mobile electron (see Fig. [T]). 
The Hamiltonian of the model looks as follows: 
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FIG. 1. A part of the spin-electron system on tetrahedral 
chain. The white (black) circles indicate the lattice sites avail¬ 
able for localized Ising spins (mobile electrons). 

In above, N is the total number of nodal lattice sites, 
'^L a Cki,a (o: =t)4.) represent fermionic creation and 
annihilation operators, respectively, while = ±1/2 
stands for the Ising spins. The parameter t > 0 takes 
into account the kinetic energy of a single mobile elec¬ 
tron delocalized over a triangular plaquette and J de¬ 
notes the Ising exchange interaction between electrons 
and their nearest Ising neighbours. Finally, Hi and He 
are the Zeeman’s terms acting on the localized Ising spins 
and mobile electrons, respectively. 

The investigated spin-electron model is exactly solv¬ 
able within the framework of a generalized decoration- 
iteration mapping transformation [^. By the use of this 
rigorous procedure, one directly derives the mapping re¬ 
lation between the partition function of the investigated 
spin-electron tetrahedral chain and the partition function 
of the uniform spin-1/2 Ising linear chain in the effective 
external magnetic field, see Eq. (2) in Ref. Q. Subse¬ 
quently, one can readily calculate the Gibbs free energy 
of the system and all important physical quantities uti¬ 
lizing the standard thermodynamic relations. 

HI. THERMODYNAMICS 

In this article, we turn our attention to the ther¬ 
modynamic properties of the spin-electron tetrahedral 
chain. Before doing so, howeve^ let us briefly men¬ 
tion the results obtained in Ref. Q. By assuming equal 
magnetic fields acting on the Ising spins and electrons 
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H = Hi = He, the model exhibits the ferromagnetic 
(FM) ground state with a full alignment of the Ising spins 
and mobile electrons into the H direction, if J < 0, or it 
passes from two-fold degenerate antiferromagnetic (AF) 
ground state with the antiparallel alignment between mo¬ 
bile electrons and their nearest Ising neighbours to the 
FM one at the field H/ J = 1, ii J > 0. In addition, 
an enhanced MCE can also be observed during the adia¬ 
batic demagnetization. More specifically, the model with 
J < 0 exhibits an enhanced MCE in the limit H/\J\ —>■ 0 
when the entropy is kept very small, while the model 
with J > 0 exhibits an enhanced MCE when the entropy 
is sufficiently close to the value S/2N = ln[(l -|- •\/5)/2]. 

Now, let us proceed to the examination of the thermo¬ 
dynamics of the spin-electron model with the ferromag¬ 
netic Ising interaction J < 0. For this purpose. Fig. [2] 
illustrates thermal variations of the specific heat for this 
particular version of the model by choosing the magnetic 
field H/\J\ = 0.001 and different values of the hopping 
term t/\J\. As one can see, C(T) curves exhibit at least a 
round high-temperature maximum, which diminishes and 
shifts towards higher temperatures upon strengthening 
t. For t/\J\ > 0.5, there also appears an additional low- 
temperature peak, which becomes more pronounced, the 
stronger the hopping term is. The observed double-peak 
structure of C'(T) curves originates from thermal excita¬ 
tions between the fully polarized ground state and the 
low-lying excited state with magnetic particles aligned 
in the opposite direction of H, which is energetically 
close to the ground state for weak fields. Indeed, the 
low-temperature maximum changes neither in the posi¬ 
tion nor height, if one considers the fixed H and the 
parameter t varies. However, when H strengthens, then 
the low-temperature peak gradually shifts towards higher 
temperatures until it merges with the high-temperature 
one (it is not shown). 
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FIG. 2. Thermal dependencies of the specific heat for the 
model with the ferromagnetic Ising interaction J < 0 at the 
fixed magnetic field H/\J\ = 0.001 and various values of the 
hopping term t/\J\. 

Another quantities, which are important for under¬ 
standing the thermodynamics, are the magnetic suscep¬ 
tibility and the total magnetization. Thermal variations 


of the susceptibility times temperature (yT) are depicted 
in Fig. [3Ka) by assuming various values of the magnetic 
field. Note that the plotted yT curves correspond to the 
particular case J < 0 and any value of t. Obviously, the 
zero-field yT product exhibits an exponential increase 
upon lowering temperature with low-temperature diver¬ 
gence that is typical for ferromagnets. By contrast, as 
far as the weak fields are considered, yT plots start from 
zero value and manifest a maximum, which flattens and 
shifts to higher temperatures with the increasing H. The 
appearance of the peak in yT curves relates to a thermal 
instability of the magnetically ordered system, whereas 
already a small temperature change necessitates a huge 
decrease of the total magnetization at weak fields, see 

Fig.ETb). 



FIG. 3. Thermal dependencies of the magnetic susceptibil¬ 
ity times temperature [Fig. EKa)] and the total magnetization 
[Fig. EJb)) for the model with the ferromagnetic Ising inter¬ 
action J < 0 at various values of the magnetic field H/\J\. 

To explore the thermodynamics of another particular 
case of the spin-electron tetrahedral chain with the an- 
tiferromagnetic Ising interaction J > 0, the temperature 
dependencies of the specific heat, magnetic susceptibility 
and total magnetization are displayed in Figs. |4] and [5] 
for different values of the hopping parameter and mag¬ 
netic field. Fig. |4] illustrates thermal variations of the 
specific heat when magnetic fields are chosen close to 
the transition value H/ J = I, where the system under¬ 
goes a zero-temperature phase transition AF-FM. As one 
sees, C'(T) curves manifest a notable double-peak struc¬ 
ture due to strong thermal excitations between ground 
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state configuration and the one of the respective low-lying 
excited state, namely from AF towards FM configura¬ 
tion (for H/J = 0.95) and vice versa (for H/J = 1.05). 
The observed double-peak structure becomes more pro¬ 
nounced, the closer the magnetic field is selected to the 
value H/J = 1 (it is not shown). 



T/J 

FIG. 4. Thermal dependencies of the specific heat for the 
model with the antiferromagnetic Ising interaction J > 0 at 
the fixed magnetic fields H/J — 0.95 and 1.05 by assuming 
various values of the hopping term t/ J. 

The aforementioned statements can be convincingly 
evidenced by thermal dependencies of the magnetic sus¬ 
ceptibility multiplied by the temperature and total mag¬ 
netization, which are displayed in Fig. [5] for various val¬ 
ues of H. Similarly as for the model with J < 0, plotted 
curves correspond to any value of t. As can be seen form 
Fig. [IKa), the quantity exponentially tends to zero 
value with decreasing temperature regardless the mag¬ 
netic field is chosen so that AF or FM phases constitute 
the ground state. The closer value of the external field 
to the AF-FM ground-state boundary is considered, the 
more pronounced temperature variation of the product 
xT can be found due to the strengthening thermal ex¬ 
citations between energetically close AF and FM con¬ 
figurations. At the critical field H/J = 1, there the 
magnetic susceptibility times temperature tends to the 
finite value yT = 0.17889 in the zero-temperature limit. 
In agreement with these observations, the total magne¬ 
tization of the system exhibits the vigorous thermally- 
induced increase (decrease) from the zero value (satu¬ 
rated value) when the applied magnetic field is chosen 
slightly below (above) the critical value H/J = 1, see 
Fig. [5Kb). The magnetization curve corresponding to 
H/ J =1 slightly decreases from the non-saturated value 
m = v^/lO ~ 0.22361 with increasing temperature. 


IV. CONCLUSION 

In this paper, we have studied the thermodynamic 
properties of the exactly solvable spin-electron tetrahe¬ 
dral chain with Ising spins localized at nodal lattice sites 
and one mobile electron delocalized over the triangular 




FIG. 5. Thermal dependencies of the magnetic susceptibil¬ 
ity times temperature [Fig. (SKa)] and the total magnetization 
[Fig. (SKb)] for the model with the antiferromagnetic Ising in¬ 
teraction J > 0 at various values of the magnetic field H/J. 

plaquette. It has been demonstrated that the double¬ 
peak structure of the specific heat curves as well as pro¬ 
nounced temperature variations of the magnetic suscepti¬ 
bility multiplied by temperature and total magnetization 
observed around the fields H/\J\ = 0 (for J < 0) and 
H/ J = 1 (for J > 0) relate to strong thermal excitations 
between ground state configuration and the configuration 
corresponding to the respective low-lying excited state. 
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